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Abstract 

The probability of a projectile nucleon to traverse a target nucleus without interac- 
tion is calculated for central Si-Cu and Si-Pb collisions. Special attention is given to 
the impact parameter range which contributes to events with large transverse energy. A 
fit to the data from E814 requires an effective NN cross section of <r e // = 54.2 ± 5.0 
mb (compared to the free space value a^ N — 30 mb ) and is interpreted as one related 
to wounded target nucleons. 
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The dynamics of a heavy ion collision is a complicated many-body problem for various 
reasons. It is the task of appropriately designed experiments to isolate one particular aspect 
of the dynamics and elucidate its physics. Wounded nucleons are one of the open problems. 
In a heavy ion collision a nucleon may undergo a sequence of collisions, which follow so 
rapidly one after another, that the nucleon is no more in its ground state, even not necessar- 
ily in any definite excited baryonic resonance ( like A or N* ). We will speak of a wounded 
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nucleon. In a next encounter with another nucleon this wounded nucleon will not interact 
with the free space NN cross section a^ n N but with an effective one a e f t . Is it possible to 
determine a e f t from experiment? 

The E814 collaboration has designed an experiment to answer this question[[j]]. At 
the energy of 14.6 GeV/nucleon, the projectile 28 Si collides with Al, Cu and Pb and the 
beam rapidity nucleons are studied as a function of the centrality of the reaction (controlled 
by a measurement of the transverse energy ). The beam-rapidity nucleons belong to the 
projectile and have not lost any energy in the reaction. Of course, in a peripheral reaction 
one always has the so called spectator nucleons, which pass the target nucleus without 
interaction. They are not interesting for our purpose. However , in a central event, e.g., 
for Si on Pb, one still sees beam rapidity nucleons. The target nucleus is transparent for 
these projectile nucleons. We expect the transparency of a heavy nucleus like Pb to be 
small. Indeed, the "survival probability" S for a projectile nucleon to pass through the 
target without any inelastic interaction has been measured to S exp — 3.5 • 10 -3 for central 
Si-Pb collisions^. Does this result contain information about wounded nucleons and 
effective NN cross sections? 

In a heavy ion collision those projectile nucleons which arrive first at the target nucleus, 
interact with ground state nucleons of the target, but may transform them into wounded 
ones. The next wave of projectile nucleons already finds wounded target nucleons and 
interacts with them via a modified, probably larger cross section and transforms them to 
a higher degree of woundedness. This procedure repeats itself for the third and follow- 
ing waves. A complicated situation like this may best be simulated by a cascade code. 
However, present day computer codes are not yet in a position to handle off-energy-shell 
situations. Instead they use on-shell baryons ( in their ground and excited states ) together 
with mesons as their basic degrees of freedom. A code like this, ARC, has been applied to 
the data from E814 by Schlagel et al. and has lead to agreement with the data without 
wounded nucleons and effective cross section, though with large error bars because of the 
low statistics. As we will argue on the basis of time scales, the basic assumptions of the 
cascade code may not be fulfilled and another approach may be better justified. It will be 
presented in this paper. 

The survival probability S for a projectile nucleon to pass through the target without 
inelastic interaction is calculated in a Glauber type approximation (straight lines, frozen 
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nucleons) as 



S(b,a eff )= Jd 2 sT p (b-s)e^ff A ^, (1) 

where b is the impact parameter of the nucleus-nucleus collision, T p (b) and T t (b) are the 
thickness functions (T(b) — J dzp(b, z), J d 3 xp(x) = 1) of projectile and target, respec- 
tively. The straight line geometry is certainly justified for the through-going nucleons. The 
use of "frozen" nucleons and the neglect of any other degrees of freedom, like mesons, 
needs justification. We consider a target nucleon and estimate the time At in its rest system 
which has elapsed between the arrival times of the first and the last projectile nucleons: 
At < 2R p /jp where 2R p ~ 7/m is the diameter of the Si projectile and 7 P ~ 15.6 is the 
Lorenz factor for the experiment under consideration. The time At « 0.5/m/c is rather 
short: (i) The target nucleon has not moved significantly in space ( frozen approximation 
is good ). (ii) According to the uncertainty principle the intrinsic excitation energy is a un- 
certain by AE > (At) -1 k, 0.4GeV and does not allow the definition of a definite excited 
state, (iii) Secondary hadrons are not yet formed since typical formation times are of order 
1 fm/c. 

The nature of a wounded nucleon is not clear. We parametrise any modification into 
an effective cross section a e ff between a ground-state nucleon (of the projectile) and a 
wounded nucleon (of the target), and determine it from experiment. Using Saxon- Woods 
parametrisations for the densities p t and p p with the surface thickness a — 0.52 fm for all 
nuclei and the half-density radius ta such that the root-mean-square radius of the nucleus 
equals the charge radius[|J, the survival probability S{b, cr e ff) is calculated for two values 
a e f f = 30 and 50mb, and the results are displayed in Fig. 1 . The experimental point, 
measured at transverse energy E£ — 15.5 GeV is also shown in the figure. It is obtained 
from the measured mean multiplicity < AI C > of beam rapidity protons by 

S(E t ) =< M c > (E t )/Z, (2) 

where Z — 14 is the number of protons in Si. We have assumed - as the authors of the 
experiment do - that the high value Ef corresponds to a central collision which is assigned 
a value 6 = 0. Then the experimental value is close to the curve a e f / = 30m6. In fact 
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the equation S(b — 0, <J e ff) — S exp (E^ = 15.5GeV) leads to a value a e ff = 31.1 ± 0.7 
mb. This result reproduces a similar calculation using uniform density distributions by the 
E814 collaboration [[!]], who have concluded that a e f j = a^ n N within error bars and no 
anomaly being visible. 

The crucial step in the argument is the assignment of b = to the central value of 
transverse energy E t . Indeed , a given value of E t determines an impact parameter b m (E t ) 
only within a band Ab(E t ). And Ab(E t ) may be very large! We study the relation be- 
tween Et, b m and Ab in the form of a probability distribution P(E t , b), which gives the 
probability that values of b contribute to events with a given value of E t . We normalize 
it as J dE t P(E t , b) = 1. With this function and the differential inelastic heavy ion cross 
section dai n /d 2 b, one can obtain the dependence of the survival function S(E t , <J e ff) as 
function of transverse energy 



qnP x Jd 2 bS(b,a eff )P(E u b)da m /d 2 b 

S{EU ^ = d^JdE t ' (3) 



d(7 in f j2 j 7-, / 77, 7 \ do"i n 



where 



dEt -j^P(E tM ^. (4) 

We will use Eq.(Q) to determine P(E t , b) from a comparison with the measured transverse 
energy distribution dain/dEt- Only then S(Et, (Jeff) can be calculated from Eqs.(Jl|) and 
(H) without ambiguity. We discuss our parametrisations. 

The inelastic cross section doi n jd 2 b is taken from the folding model 

(kTin {b) = 1 - exp[-ag N A p A t J d 2 s T p (b- s)T t (s)]. (5) 



d 2 b 

For central collisions (small value of b ), the exponential is practically zero and any uncer- 
tainties, e.g., the choice of a^ n N , are unimportant. The parametrisation of P(E t , b) is more 
model dependent. We choose a Gaussian parametrisation 



P(E t ,b) = s^^ expi- 1 \ *V J }, (6) 



i r [E t -E t {b )f 

^m eXH 2of(6) 

which satisfies the normalization condition. We make the usual assumptions for the 
functions E t (b) and at(b) 
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E t (b) = N(b)e , 



(7) 



(8) 



Here N(b) is calculated either in the " collision model" A^ C (6)|Q] 
cleon model" iV^, (&)[§], 



or in the "wounded nu- 



N c {b) 



= o, 




(9) 



N w (b) = A t / d 2 bT p (s)exp[-ag N T t (b - s)} + (t «— > p). 



(10) 



N c (b) equals the mean number of NN collisions in a projectile-target interaction with im- 
pact parameter b, and N w (b) equals the number of nucleons in the overlap volume of pro- 
jectile and target. The proportionality constant eq between the observed transverse energy 
E t and N(b) depends on the dynamics of hadron production, but also on the experimental 
set up (chosen rapidity interval and acceptance of counter). It will be a fit parameter. In 
the collision model, which we will use, only the product eocr^ enters and the physics of 
Eq.(||) is completely independent of the choice of af n N . For the wounded nucleon model 
this is only approximately true. The ansatz of oc N(b) corresponds to the hypothesis that 
we deal with statistical fluctuations in N(b). It is known, however, that the proportionality 
constant ui in Eq.(||) depends strongly on the rapidity interval of the accepted particles [Q]. 
The reason is not clear[^]. We take u> as a free parameter. 

Using expressions (||) to (|8]) we have calculated d<Ji n /dE t and have varied the pa- 
rameters eo an d tJ until the data for Si-Pb are fitted. The result is shown in Fig. 2 with 
ag N = 30mb, e = 0M7GeV, uj = 7.8 for the collision model and e Q = OAGeV, oj = 6. 
for the wounded nucleon model. Both models for the calculation of N(b), describe the 
data Si-Pb equally well. But when the same parameters are used to describe dtJi n /dE t for 
Si-Al and Si-Cu collisions, the collision model gives better fits and therefore we discard 
the other model. 

According to Eq.(Q), P(E t , b) determines the integration region in impact parameter 
b, which contributes to the integral for a given value of E t . As an example, Fig.l shows 
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the weight function P(E%, b). Its width is unexpectedly large. In particular the survival 
probability S(b, <r e ff) changes considerably in the 6-region, where P(E^, b) is large. This 
invalidates the assumption to put the experimental point at b = 0. 

Since P(Et, b) plays such important role in the understanding of the experiment, we 
have studied its properties in more detail. For the functions P(E t ,b) used to fit dai n /dE t 
we have calculated the value b m (E t ) of the impact parameter which contributes maximally 
to a given E t and Ab(E t ) which is the width of the band of impact parameters around 
b m {Et)- We define these two quantities by 



db 



P(E t ,b)\ b = bm(Et) =0, (11) 



and 



[Ab(E t )} 2 = J d 2 b(b~b m ) 2 P(E t ,b)/ J d 2 bP(E u b). (12) 
For the Gaussian parametrisation Eq.(||) one has approximately the relations 

E t [b m {E t )] = E u (13) 



where N'(b) = dN/db. Note that the width of the band in impact parameter , which 
contributes to a given E t is proportional to the width oj of the Gaussian of P(E t , b). 

Fig. 3 shows the position b m (E t ) of the maximum of P(E t , b) and the width Ab(E t ) 
as a function of E t for Si-Pb collisions. As expected the position b m goes to zero and 
then stays zero for large values of E t . However, the width Ab shows a break at E® where 
b m (Et) becomes zero. This is related to a sudden change in the shape of P(E t , b) at E®. 
For E t > E® the width Ab decreases only very slowly. 

Fig.l shows as example the distribution P(E t , b) for a central collision ( £" t c = 15.5 
GeV ). We have compared this form with the corresponding one Pu(Ef, b) calculated [j|] 
with the help of the code UrQM D\J^. The shape of Pjj is similar to P, but narrower. 
(The value of < b > is smaller by about 25%). However, also da/dE t calculated from 
UrQMD falls off faster than the experimental data. Using the functions S(b, <r e ff) and 
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P(Ef, b) as shown in Fig.l, the calculated value of S(E%, cr e ff) gets contributions from a 
considerable range of values of b. We find 

S{El a eff )/S(b = 0, a eff ) = 2.86 (15) 

for a e ff = 30mb. Since S(b — 0, a e ff = 30m&) corresponded essentially to the experi- 
mental value, S(E£, a e ff — 30m6) does not. A solution of the equation for <7 e // 

S(E c t ,<j eff ) = S exp {E r t ) (16) 

leads to 

a eff = 54.2 ± 5.0 ± 6.4 ± 5.6m6. 

The meanings of the errors are explained below and in the caption of Table 1 . The same 
equation for central events ( E^ = 9.5 GeV )in Si-Cu collisions leads to a e jf = 50.3 ± 
15.0±4.6±6.4m& , in agreement with the value from Si-Pb. Since Al is a smaller nucleus 
than Si, the beam rapidity nucleons from Si-Al will always be contaminated by spectator 
nucleons and not very sensitive to effects of the nuclear transparency. The results of the 
values for a e ff are compiled in Table 1. We have also calculated a e ff by using Pu(E^r, b) 
from U rQMD and have obtained a e f / — 40.3 ± 3.0m6, where the error reflects only the 
statistical fluctuations in Pjj. 

The result for a e f f from our analysis is at variance with the conclusion by the authors 
of the E814 experiment [[IJ and of the cascade calculation of Ref. [Q]. While we are unable 
to understand the difference to the cascade calculation^, the difference to the argument 
of the E814 collaboration is clear: They assume that central collisions means 6 = with 
very small width, while our analysis shows that even for central collisions a fairly wide 
band of impact parameters contributes. The width of this band depends on the parameters 
eo and uj in the function P(E t , b). How sensitive are the data to these parameters? Table 1 
gives the uncertainties in the extracted values of a e tf if one assumes uncertainties of 5% 
in eo and 10% in ui, respectively. A reduction by 25% in uj would reproduce the results of 
UrQMD for da/dE t and a e ff but would be in contradiction to the data for da/dE t and 
therefore the extracted value <j e / / ~ 40m6 is doubtful. According to the results presented 
in Table 1 and from the cascade calculations the extracted effective cross section seems to 
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be definitely larger than the inelastic cross section af n for a NN collision in vacuum. 

The difference may be interpreted in terms of wounded nucleons. We try to estimate the 
effect by the following model [|l(J which a wounded nucleon N + consists of a nucleon plus 
a semi-hard prompt gluon, which is radiated as bremsstrahlung in quark-quark interactions 
during an inelastic NN collision JTT|]. We furthermore use the additive quark model ( a 
nucleon consists of 3 quarks ) to relate the NN cross section to the nucleon-quark (Nq) 
one via a NN — 3a Nq . Then 

If in addition one assumes a Ng /a Nq rs 9/4, the color factor, one obtains cr^ n N+ w 50mb. 
In view of all the uncertainties in the definition of <7 e //, in the model for N + and in the 
additive quark model, the close agreement between experiment and the theoretical estimate 
must be considered fortuitous, but it may show that the order of magnitude of the extracted 
value for the effective cross section is not unreasonable. 
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Captions: 

Fig.l Survival probabilities S(b, <J e ff) of a beam rapidity nucleon calculated for a e ff — 
30 mb and a e // = 50 mb (solid lines) and the correlation P(E^ = 15.5GeV, b) between 
transverse energy and impact parameter (dashed line) for Si-Pb collisions as a function of 
b. The data point is obtained from the measured multiplicity < M c > = 0.049 ± 0.005 of 
beam rapidity protons at E$ . 

Fig.2 Transverse energy distributions in the rapidity range —0.5 < r\ < 0.8 for a Si 
beam at 14.6GeV/nucleon on different targets. Data are taken from [7]. The calculations 
refer to the collision model (solid lines) and the wounded nucleon model (dashed curves). 
For each model one set of parameters e , u) has been determined by a fit to the Pb data and 
then applied to the other data. 

Fig.3 The position b m {E t ) of the maximum of P(E t ,b) in impact parameter b (solid 
line ) and the width Ab(E t ) around it (dashed line ) as a function of E t for Si-Pb collisions. 

Table 1 The fitted effective NN cross sections from the experimental value of the mul- 
tiplicity < M c > of beam rapidity nucleons in collisions of Si-Pb and Si-Cu at AGS 
energy of 14.6GeV/nucleon. If one assumes that the El corresponds to b = and uses dis- 
tributions of uniform density (u.d.) or Saxon- Woods (S.W.) one finds the values in columns 
4 [1] and 5, respectively. The last column corresponds to our result. The first error cor- 
responds to the uncertainty in < M c >, while the second and third errors correspond to 
uncertainties of 5% in e an d 10% in ui, respectively. 
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££(GeV) 


< M c > 


cr e// (mb) 


b=0(u.d.) b=0(S.W.) E t = E c t 


Pb 


15.5 


0.049±0.005 


28.8±0.5 31.1±0.7 54.2±5.0±6.4±5.6 


Cu 


9.5 


0.66±0.09 


28.8±1.8 35.8±2.5 50.3±15.0±4.6±6.4 



Table 1 
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